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ABSTRACT 






In an effort to Improve the current solutions In the design and analysis of 
liquid propulsive engines, a computational fluid dynamics (CFD) model capable of 
calculating the reacting flows from the combustion chamber, through the nozzle to the 
external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea 
level, was Investigated as a sample case. The CFD model, FDNS , Is a pressure based, 
non-staggered grid, vlscous/lnvlscld. Ideal gas/real gas, reactive code. An adaptive 
upwlndlng differencing scheme Is employed for the spatial discretization. The upwind 
scheme Is based on fourth order central differencing with fourth order damping for 
smooth regions, and second order central differencing with second order damping for 
shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a 
PARASOL finite rate chemistry algorithm using the point implicit method. In this 
study, the computed flow results and performance compared well with those of other 
standard codes and engine hot fire test data. In addition, the transient nozzle 
flowfleld calculation was also performed to demonstrate the ability of FDNS In 
capturing the flow separation during the startup process. 

INTRODUCTION 

Several limitations have been found with the current industrial standards in 
calculating the flowfleld and performance of liquid rocket engines. For example, the 
Method of Characteristics (MOC) was generally used to analyze the supersonic nozzle 
flowfleld; however, factors that could influence the numerical results included 
transonic flow analysis, the location and flow properties of the start line, along 
with the number and distribution of points on the start line. Moreover, MOC is 
restricted to supersonic flow. Composite codes and decoupled solutions have to be 
resorted to resolve the wall boundary layer, the combustion chamber flow, and the 
exhaust plume. Limitations in axisymmetrlc formulations and steady-state solutions 
are two other shortcomings. The recent advancement of CFD technology has shown 
potentials in transient flow simulations (ref. 1), flow reversal prediction 
capabilities, and a unified solution (ref. 2) for subsonic, transonic and supersonic 
flowfields. Furthermore, general two dimensional and three dimensional geometry can 
be described; and applications to dual throat engines, scarfed and Integrated 
nozzles, and cluster nozzles may be achieved. This study produced a CFD model for 
the description of the unified liquid rocket engine flowfields, and the prediction of 
the thrust performance with the unified solution. Code validation was achieved by 
choosing SSME firing at 100* power at sea level as the sample case, and by verifying 
the code from such options as ideal gas, invlscid, non-reacting, adiabatic wall to a 
more complicated real gas, viscous, equilibrium, and specified wall temperature 
calculation. The computed results were compared to those of various industrial CFD 
codes and hot fire experiments. Finally, a transient calculation was performed to 
demonstrate the CFD code's ability in capturing the shock evolution and the flow 
separation during the transient operation. 

GOVERNING EQUATIONS 

The basic equations employed to describe a unified liquid rocket engine 
flowfleld are the axisymmetrlc, multi -component transport equations. The 
classification of the governing equations changes from one point to another in the 
unified flowfleld. That is, elliptic for subsonic flow, parabolic for sonic flow, and 
hyperbolic for supersonic flow. A generalized form of these equations written in 
curvilinear coordinates is given by 

(i/j)(3pq/at) = at-pUfl + i/G 1J (aq/ae j )]/a< 1 + s q (D 

where J, U,* and G^j are the Jacobian of coordinate transformation, the contravariant 
velocity and the geometrical matrices, q represents 1, u, v, h, k, e and p^. 
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respectively. These are equations of continuity, x and y momentum, enthalpy, 
turbulent kinetic energy, turbulent kinetic energy dissipation rate, and density 
fractions. The effective eddy viscosity u is determined by solving a extended two- 
equation k-e turbulence model (refs. 3 and 4). Source terms S q are given by 


S q = l/J 


0 

-Px + v [ v (u j) x ] - 2/3 { wu j) x 
-Py + ^t^(Uj)y] - 2/3 ( WUj) y 
Dp/Dt + ® + £ J rp p n ?T - Ih^ 

p(P r - e) 

p(e/k) (C^p-Cje+C^^/e) 
n = 1 f t • • f N 


( 2 ) 


where 0 is the energy dissipation function, J n is diffusion flux, C pn is heat 
capacity, and W n denotes mass production rate for species n. P r represents the 
turbulent kinetic energy production term. C 1# C 2 and C 3 are model constants for the 
two-equation turbulence model. 

NUMERICAL SCHEMES 

An adaptive upwinding scheme was employed to approximate the convective terms of 
the momentum, energy and continuity equations; the scheme is based on second and 
fourth order central differencing with artificial dissipation. First order upwinding 
is used for species and turbulence equations, since the parameters involved are 
positive quantities. Different eigenvalues are used for weighing the dissipation 
terms depending on the conserved quantity being evaluated, in order to give correct 
diffusion fluxes near wall boundaries. Adding the dissipation term to the convective 
fluxes F in computational coordinate £ produces 

3F/ae - (F i+1 - F ■}-,) / 2 - (d 1+1/2 - d M/2 ) (3) 

The dissipation term, d, is constructed such that a fourth-order central and 
fourth-order damping scheme is activated in smooth regions, and a second-order 
central and second-order damping scheme is used near shock waves. Since the Jacobian 
matrices of the Euler fluxes have eigenvalues of U, U+c and U-c, it is sufficient to 
use the magnitudes of these eigenvalues to weigh the dissipation terms to maintain 
the smoothness of the solution without losing accuracy. |U|+c was used for the 
continuity equation and |U| was used for other transport equations in FDNS. A 
general form of the dissipation term is given by 

dj+ 1/2 = 0.25[2ej|pU| + e 2 (|U| + c) ] i+i/ 2 < <3 i+i “ 

+ [ e 3 ( 1 - e i)MAX[ 0 . 5f p ( | u | + |v|),2|pU|] 

+ e 4 (|U| + c) ] i+i /2 ( q 1-1 - 3qi + 3q 1+1 - q i+2 ) (4) 

where a i = |p 1+1 - 2p^ + p H l/(Pi+i + + Pi-l) 

^ i+1/2 = MAX(a^, a 1+ i) 

dj = MAX [ A , MIN (1 . 0 , 2 5 0 i + i/ 2 ) ] 

d 4 « MAX ( 0.0, 0.01 - O.250 i+1/2 ) 

Different values for dissipation parameters e j, e 2 , e 3 and e 4 are used for the 
continuity, energy and momentum transport equations. f represents a local cell flow 
area. Unity X corresponds to a fully upwinded differencing scheme and vanishing X 
corresponds to a central differencing scheme for the convective terms in smooth 
regions. Table I summarizes these parameters. 
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Table I . Dissipation Parameters 



Momentum & Energy 

Continuity 

e 1 


0 

e 2 

0 

9 i+t/2 

«3 

0.015 

0 

e 4 

0 

<*« 


A pressure based method was selected such that a wide range of flow speeds could 
be analyzed with the same code. Successful results of viscous flow computations 
using pressure based methods have been reported (refs. 5 and 6). For high speed flow 
cases, a hyperbolic pressure correction equation was employed by perturbing the 
density in the mass conservation equation. This provides a smooth transition from 
low to high speed flow characteristics. For time accuracy, a time-centered, time- 
marching scheme with a multiple pressure corrector algorithm was employed. In 
general, a noniterative time-marching scheme was used for time dependent flow 
computations (ref. 1); however, subiterations can be used if necessary. The multi- 
corrector procedure will be described below. 

A simplified momentum equation was combined with the continuity equation to form 
a pressure correction equation. The simplified momentum equation can be written as: 

apu/at = - vp* 
or, in discrete form, 

u' 1 = - 0( At/p) vp* (5) 

where 0 represent a pressure relaxation parameter. The velocity and density fields 
in the continuity equation are then perturbed to form a correction equation. That 
is, 

[ap/vt + v(pu 1 )] ft+1 = a(p n +p')/9t + (p n +p» ) (u n 1 +u f i ) ] = o 

By neglecting the p'u* terms, the following equation results. 

ap'/at + v(u 1 p l ) + v(pu , i ) = - [ap/at + vipu^)] 0 (6) 

Substituting Eq. (5) into Eq. (6) and letting p 1 = p'/RT, the following pressure 
correction equation is obtained. 

a(p'/RT)/at + 7 [(u/rt)p'] - v(0At vp 1 ) = - [ap/at + vfpuT)]" (7) 

To provide smooth shock solutions the adaptive dissipation terms described above were 
added to the right hand side of Eq. (7). Once Eq. (7) is satisfied, the velocity 
field and the pressure field are updated through Eq. (5) and the following relation. 


The density field is then updated by applying the equation of state. To ensure that 
the updated velocity, density and pressure fields satisfy the continuity equation, 
the above pressure correction solution procedure is repeated several times before 
marching to the next time step. This represents a multi-corrector solution 
procedure . 

The chemistry source terms were evaluated with a point implicit procedure before 
the species equations were solved. The determination of equilibrium chemistry source 
terms were based on the CHMQGM algorithm (ref. 7). 

CHMQGM EQUILIBRIUM CHEMISTRY ALGORITHM 

The chemistry source terms must be computed for each grid point at each time 
step. The point implicit procedure allows the equilibrium calculation for each grid 
point to be independent from that of its neighbors and may be treated in isolation. 
Since kinetics loss was estimated to be negligible in the SSME thrust chamber and 
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nozzle (ref. 8 ) f equilibrium Is assumed to exist for the entire computational domain, 
including the exhaust plume region. 

In general, if there are N distinct chemical species composed of M chemical 
elements, then the algebraic system to be solved consists of N-M nonlinear 
equilibrium equations and M linear element -conservation relations. Let A^ be 
the chemical symbol of the ith species. While a 1s and b is stand for the 
stoichiometric coefficients for the forward and backward reactions. The equilibrium 
reactions are of the form 

N N 

I a 1s Ai = £ b is A,, s=l , . . .N-M 

1=1 1=1 

The equilibrium reactions lead to algebraic relations of the form 
N N 

K s n [c t ] is o n [c^is, s=l , . . . N-M ( 8 ) 

1=1 1=1 

where K s is equilibrium constant for sth reaction. The linear element-conservation 
reactions are of the form 

N 

2 d 1jCl = Cj j=l,...M (9) 

1=1 

where Cj is species molar concentration, Cj denotes element molar concentration and 
d^j represents the number of element in a species. 

CHMQGM algorithm has shown that a reduced system can be obtained by substituting 
the M linear relations, Eq.(9), into the N-M nonlinear relations, Eq.( 8 ). As a 
result, a iterative Newton-Raphson technique was used to solve a system of N-2M 
equations rather than the larger equivalent NxN system. The equilibrium constant 
approach of CHMQGM is not as general as a minimization of Gibb's free energy approach 
(ref. 9) for an arbitrary chosen chemical system. However, for a generally well- 
known hydrogen/oxygen system such as the one being used in SSME, CHMQGM appears to be 
more efficient. 


HYDROGEN/OXYGEN EQUILIBRIUM CHEMISTRY 

The Hydrogen/Oxygen equilibrium reactions used in this study is a subset reduced 
from the hydrocarbon combustion partial equilibrium system (ref. 7). It consists of 
six species and four reactions as shown in Table II. 


Table II. Species and reactions considered 


Subscript 

Species 

Subscript 


Species 


1 H20 

Water vapor 

4 

OH 

Hydroxyl 

radical 

2 H 

Atomic hydrogen 5 

H 2 

Molecular 

hydrogen 

3 0 

Atomic oxygen 

6 

0 2 

Molecular 

oxygen 

O 

CM 

II 

CM 

O 

K|C 6 = 

(c 3 ) 2 


( 10 ) 


X 

CM 

II 

rs» 

X 

K ^ 5 = 

(C 2 ) 2 


( 11 ) 


20H = 0 2 + H 2 K 3 (c 4 ) 2 = 

C(jC 5 


( 12 ) 


2H fi = 2H 2 + 

0 2 K,(c ,) 2 = 

(c 5 ) 2 c 6 


(13) 



The element-conservation equations are 
C 0 — c ] + c 3 + c 4 + 2 c 6 
C# = 2 c j + cj + C 4 + 2 c 5 

Substitute Eq.(10) and Eq.(ll) to Eq.(12) and Eq.(13), we have 


c 6 = Ri(c 3 ) 4 




1/K t 


( 14 ) 

(15) 


(16) 
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Rj(c 2 ) 2 

R 2 - 

l/Kj 

R 3 C 2 C 3 

*3 - 


RjCcpj 

R 4 - 

1 / ( K £ 4 ) 1/2 


(17) 

(18) 
(19) 


Substitute the linear term in equations (16-19) into the linear element- 
conservation equations (14 and 15), yielding two nonlinear equations and three 
unknowns in C 2 » c 3 , and c§. 


Co = 


R<c 5 C 3 + C3 + + 2 H 1 (c 3 ) 


2 


( 20 ) 


C H 


2R4C5C3 + c 3 + R ^^3 + 2R2(c 2 ) 


2 


(21) 


The final reduced equations can then be obtained by substituting Eq.(20) to 
Eq . (21) and to eliminate c 5 in Bq.(20) with Eq.(17), we have 

R^ 4 (c 2 ) c 3 + R 3 PJC 3 + 2Rj(c 3 ) 2 + c 3 - C 0 = 0 


and 


R3 c 2 c 3 **" 4R^(c 3 ) + 2 c 3 — 2R 2 ( C 2 ) “ C 2 + (C^-2Cq) 0 0 

This reduced system consists of one cubic equation and one quadratic equation 
with two unknowns, and can readily be solved by Newton-Raphson 1 s iterative method. 

BOUNDARY CONDITIONS 

Fixed chamber conditions were used at the inlet of the combustion chamber. 
Subsonic boundary conditions were used at the inlet of the external ambient air and 
the pressure was extrapolated. This is to allow the transient disturbances sent from 
downstream to permeate through the boundary. Flow properties at the wall, centerline 
and exit were extrapolated from those of the interior domain. To obtain a unique 
solution for the SSME sample case at sea level, a fixed pressure was applied to the 
outermost point of the downstream ambient exit boundary. A wall temperature 
distribution (ref. 2 ) was specified for the viscous, real gas, and regeneratively 
cooled wall calculation. For transient thrust chamber and plume flowfield 
calculations, a digital transient model simulation ( DTM ) was used as the combustor 
upstream boundary conditions. 

SAMPLE CASES FOR SSME AT 100* POWER LEVEL 

SSME operating conditions at 100* power level and sea level were used for the 
calculations. The actual values used are shown in Table III. 


Table III. Operating Parameters 


Mixture Ratio 


6.000 

02 inlet temperature 

(K) 

90.56 

H2 inlet temperature 

(K) 

20.56 

Chamber Pressure (PSIA) 

2935.7 

Chamber Temperature 

<K) 

3639.0 

Geometric Area Ratio 


77.5 


RESULTS 

INVISCID, IDEAL GAS , AND ADIABATIC WALL CALCULATION 

FDNS calculations were made for the inviscid, ideal gas and adiabatic wall 
conditions in order to compare to the equivalent M0C solutions. The domain of 
computation was started from the injector faceplate, through the throat and ended at 
the nozzle exit plane. Computational results obtained from running the PARC code 
(ref. 10) was also used for comparison. Cases for different specific heat capacity 
ratios ( 7 = 1 . 14, 1.1875 and 1.25) were performed. Fig. 1 shows the comparison of 
vector and Mach number contour for FDNS and PARC solutions at 7=1.1875. The vectors 
showed a transition from subsonic flow in the main combustor to sonic flow at the 
throat, and to supersonic flow in the nozzle. While FDNS used a pressure based 
method, PARC used a density based procedure. Both codes have captured the nozzle 
shock stemming downstream from the throat, and exhibited similar Mach number contours 
throughout the computational domain. A sonic line can be observed near the throat 
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region. The computations were made on the same grid with a grid size of 111x65. 

Fig. 2 shows the comparison of centerline and wall pressure distributions for 
7 = 1 . 1875 . The FDNS predictions compared well with those of the MOC and PARC results. 
The sonic start lines obtained from FDNS calculations were used for MOC calculations. 
The computational time for a typical FDNS ideal gas calculation was estimated as 
1.03E-4 CPU seconds per grid per step on a NASA/MSFC CRAY-XMP. Five hundred 
iterations were required for approximate convergence. 

Performance calculations showed specific impulses (ISP) of 426.6, 452.5, and 
472.3 for 7—1 . 25 , 1.1875 and 1.14, respectively. Real gas calculations indicated 
that the heat capacity ratio ranged from 1.14 at the combustor inlet, to about 1.25 
at the nozzle exit. The thermodynamic properties of the combusting gases are very 
important parameters to the accurate prediction of the liquid rocket engine 
flowfields and performance. A real gas thermodynamics calculation is therefore in 
order . 

INVISCID , REAL GAS , AND ADIABATIC WALL CALCULATION 

The computational domain was extended to include the exhaust plume, and the 
ambient air for the low altitude FDNS real gas calculations. Seven species were 
considered. Including H 2 , 0 2 , Hp, 0, H, OH and N 2 . Chemical equilibrium was assumed 
to be valid for the entire flowfield. Fig. 3 shows the comparison of FDNS and MOC 
predicted centerline and wall Mach numbers. MOC used equilibrium properties table 
obtained from CEC thermal equilibrium calculations (ref. 9). The FDNS predicted 
centerline Mach numbers agreed well with those of the MOC prediction. The FDNS 
predicted wall Mach number matched MOC calculation from approximately the center of 
the nozzle to the exit plane, but underpredicted near the throat. The 
underprediction may have been caused by two reasons: First, a constant perpendicular 

sonic start line (Mach number=1.01) was assumed for the MOC calculation that was 
different from the FDNS calculated sonic line; second, the artificial dissipation 
used to capture shock by FDNS may have generated some total pressure loss, near the 
starting point of the nozzle shock at the wall, and as a result, the underprediction 
of the Mach number near the throat. Nevertheless, the Mach number comparison for the 
FDNS and MOC was reasonably good. Fig. 4 shows the comparison of centerline and wall 
pressure distributions. The FDNS predicted centerline and wall pressures agreed well 
with those of the MOC predictions. RAMP code (ref. 11), overpredicted the centerline 
pressure slightly near the exit plane. 

The FDNS predicted SSME nozzle specific impulse is compared to that of MOC and 
TDK predictions in Table IV. TDK calculation (ref. 8) used an eight-reaction finite 
rate kinetics mechanism. The closeness of all three ISP values indicated that the 
kinetics loss and shock-kinetics interaction is probably negligible under these 
circumstances. The computational time for a typical FDNS real gas calculation was 
estimated to be 2.36E-4 CPU seconds per grid per step. Approximately Twenty-five 
hundred iterations were required for near convergence. The grid size used in this 
calculation was 201x81. 


Table IV. ISP comparison 


Code 

ISP 

Chemistry 

FDNS 

460.4 

CHMQGM 

MOC 

458.4 

CEC 

TDK 

459.2 

Finite Rate 


VISCOUS, REAL GAS , AND COOLED WALL CALCULATION 

The simulation of 100* power level SSME fired at sea level, was completed by 
adding the viscous terms and a specified wall temperature to the FDNS calculation. 
The extended k-e turbulence model and wall function approach were used to calculate 
the viscous flow. The specified wall temperature was obtained from ref. 8. Fig. 5 
shows the Mach number contour of a FDNS calculated conical nozzle flow. The 
operating conditions were identical to that of the bell-shaped SSME nozzle as 
indicated in Table III, including the nozzle length and the area ratio. This was to 
see if the pressure based FDNS CFD model can capture a perfect normal shock from a 
nearly one dimensional nozzle flow. From Fig. 5, a nozzle shock formed after the 
inflection point, hit the centerline, and reflected inside the nozzle; a lip shock 
formed at the nozzle exit and extended to the triple point in the plume region, 
where a Mach disc dropped beneath the triple point and perpendicular to the axis. 
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The computed Mach number contour of the overexpanded SSME nozzle flowfleld Is 
shown In Fig. 6. Unlike the one-dimensional conical nozzle flow, a distinctive two- 
dimensional flow behavior is observed. A slightly curved Mach disc developed in the 
plume region; the chamber total temperature was recovered behind the Mach disc. The 
disc location was estimated to be 1.1 nozzle exit radii from the exit plane, with a 
width of about 1.1 nozzle exit radii. The location and width of the computed normal 
disc agreed well with those of the hot fire test (ref. 12). The thrust performance 
comparison of the experimental data, the FDNS and TDK/BLM results are shown in Table 
V. BLM is a boundary layer code that provides TDK with the boundary layer solution. 
The FDNS calculated ISP value compared very well with that of the hot fire test (ref. 
13) and the TDK/BLM calculation. 


Table V. ISP comparison 



ISP 

Chemistry 

Experiment 

453.3 


FDNS 

453.4 

CHMQGM 

TDK/BLM 

452.3 

Finite Rate 


The typical grid size used in this calculation was identical to that of the 
inviscid case. The computational time for one typical run was 2.77E-4 CPU seconds 
per grid per step. Approximately Three thousand iterations were required for near 
convergence. The CPU time for a nonreacting case was 2.32E-4 seconds per grid per 
step. That means only about 20 % increase for a equilibrium chemistry calculation. 

TRANSIENT SSME THRUST CHAMBER AND PLUME FLOWFIELD CALCULATION 

The existing SSME thrust chamber was designed to give optimum performance under 
the constraint that the nozzle lip pressure is fixed around 5.5 PSI. This pressure 
was selected to avoid flow separation at sea level during full power level. The 
vector plot of CFD calculations showed no flow separation under these operation 
conditions, as presented in Fig. 1. Asymmetric flow separation during transient 
operation is then the phenomena of prime interest, since the resultant side forces 
could be highly detrimental to the engine. As a first attempt, a transient 
calculation was performed on a axisymmetric grid to study the flow separation inside 
a nozzle. 

The transient thrust chamber and plume flowfleld was calculated as a viscous, 
real gas, and reacting (CHMQGM) flow calculation. The regeneratively cooled wall 
temperature distribution was assumed to be the same as that of the steady state case. 
The grid domain used was also identical. To simulate the actual hot fire operation, 
DTM simulation was used to generate the transient upstream boundary conditions. 

These transient upstream flow properties included temperature, pressure, and 
propellants flow rates at the injector end. The igniter flow and the coolant flow 
were combined with the main injector flow to account for the total flow rate. Sea 
level pressure and temperature were used as ambient properties. 

Pressure contours at different time cuts are shown in Fig. 7a-7f. Fig. 7a shows 
the pressure contour at 1.58 seconds. It can be seen that a normal shock wave was 
moving down the thrust chamber. This strong shock imposed a large adverse pressure 
gradient in the nozzle. As the boundary layer was not able to negotiate this adverse 
pressure gradient due to wall friction, the flow separated at the foot of the rear 
oblique shock. Another front oblique shock stemmed from the triple point, albeit it 
is somewhat smeared by the encroachment of the back flow. A nozzle shock can be seen 
developing from the inflection point near the throat. Fig. 7b shows the progress at 
2.13 seconds. The strong shock moved down some more and the flow separation 
persisted. Fig. 7c occurred at 2.53 seconds. The developing nozzle shock that hit 
the centerline may have strengthened and bent the normal shock. From this time to 
about 2.7 seconds, that is when the normal shock and the two oblique shocks were 
about to move out of the nozzle exit plane, the back flow may have formed a 
triangular formation that sometimes was referred as the "teepees". In Fig. 7d, the 
normal shock has left the nozzle and developed a Mach disc at 3.17 seconds. This is 
consistent with the observation. The lip shock and the reflected shock were in the 
making. They were further developed at 3.50 seconds, as shown in Fig. 7e. The 
leading pressure wave has since arrived and past the rightmost grid domain. 

The transient operation takes about 5 seconds to reach steady state. Fig. 7f 
shows the steady state pressure contours of the SSME firing at the sea level. This 
figure represents the same flowfleld as that of Fig. 6. There is no flow separation 
observed at this point. The total CPU time was approximately 122 hours. 
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CONCLUSIONS 


A time accurate, pressure based reactive CFD model FDNS has been developed to 
analyze a unified SSME liquid rocket engine flowfleld. The flow domain included the 
combustion chamber, nozzle, exhaust plume and ambient. The code was validated 
through several systematic stages. For example, the FDNS predicted invlscld nozzle 
flowfields were first compared to those of the MOC, PARC, and RAMP solutions. The 
comparisons showed that the FDNS predictions were reasonably accurate for invlscld 
ideal gas and real gas cases. The viscous, real gas and specified wall temperature 
case have also shown excellent comparison in terms of Mach disc width and location to 
those of the hot fire test. Furthermore, the FDNS predicted SSME ISP values were in 
excellent agreement with those of the MOC and TDK predictions, and the Rocketdyne hot 
fire test data. A ax i symmetric flow separation was also Investigated through the 
transient thrust chamber and plume flowfleld calculation. This first attempt 
demonstrated the capability of FDNS in studying the flow separation and shock 
evolution in the engine startup process. 
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FLg.l Velocity vector and tnach number contours 
for the Invlscld, Ideal gae (1-1.1075) and 
adiabatic wall SSME nozzle flowfleld 



Fig. 2 Comparison of pressure distributions for 
the Invlscld, Ideal gas (1-1.1875) and adiabatic 
wall 5 SHE nozzle flowfleld 


□ 



Fig. 3 Comparison of Mach numbers for the 
Invlscld, real gas and adiabatic wall SSME 
nozzle flowfleld 



Fig . 4 Comparison of pressure distributions for 
the Invlscld, real gas and adiabatic valL SSME 
nozzle flowfleld 
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Fig. 7 Transient SSME pressure contour evo 
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Fug. 6 The Mach number contours for the viscous 
real gas and specified wall temperature SSME 


flowfleld 





